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Voltage-flux-characteristics of asymmetric dc 

SQUIDs 

Jochen Miiller, Stefan Weiss, Rudolf Gross, Reinhold Kleiner, Dieter Koelle 



Abstract — We present a detailed analysis of voltage-flux 
V (^-characteristics for asymmetric dc SQUIDs with vari- 
ous kinds of asymmetries. For finite asymmetry a/ in the 
critical currents of the two Josephson junctions, the minima 
in the V (^-characteristics for bias currents of opposite po- 
larity are shifted along the flux axis by A<J> = aiPi, relative 
to each other; /3j, is the screening parameter. This simple 
relation allows the determination of ctj in our experiments 
on YBa2Cii307_<5 dc SQUIDs and comparison with theory. 
Extensive numerical simulations within a wide range of 
and noise parameter F reveal a systematic dependence of 
the transfer function Vj, on ctj and ctji (junction resistance 
asymmetry) . As for the symmetric dc SQUID, V$ factor- 
izes into g(T(3^) ■ f(cti,0r,), where now / also depends on ctj. 
For (3l~5 we find mostly a decrease of V$ with increasing ctj, 
which however can only partially account for the frequently 
observed discrepancy in V$ between theory and experiment 
, for high-T c dc SQUIDs. 

Keywords — High-temperature superconductors, SQUIDs, 
superconducting devices. 



I. Introduction 

THE observation of a significant discrepancy between 
numerical simulations and experimental results ob- 
tained for direct current (dc) superconducting quantum 
interference devices (SQUIDs) based on high-transition- 
temperature superconductors (HTS) is one of the most im- 
portant unsolved problems for HTS dc SQUIDs which seri- 
ously hinders their optimization for applications (T). HTS 
dc SQUIDs show frequently asymmetric behavior which 
may be attributed to the large spread in the critical current 
i Iq and normal resistance R of HTS Josephson junctions. 
This may lead to asymmetric critical current I c or voltage 
V vs. external flux $ characteristics of the dc SQUID Q, 
1 and can affect the transfer function V$ = \dV/d$\ max 
which is defined as the maximum slope of the U( ( I > )-curves. 
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However, such an asymmetry has been usually neglected in 
numerical simulations of V$ for HTS dc SQUIDs. 

In this paper we present a detailed study of the impact 
of asymmetry on the U($)-characteristics and in particu- 
lar on the transfer function of dc SQUIDs. We first intro- 
duce the main parameters which define the asymmetric dc 
SQUID (Sec.||). Then we show that an asymmetry in Iq of 
the two Josephson junctions leads to a shift of the I c (&)- 
and in the V r ($)-characteristics, which can be used to de- 
termine the critical current asymmetry experimentally, as 
demonstrated on dc SQUIDs with YBa 2 Cu 3 7 -5 (YBCO) 
bicrystal grain boundary Josephson junctions [|J (Sec. [II). 
We present numerical simulation results for V$ which we 
obtained within a wide range of parameters including the 
limit of large thermal fluctuations which are important for 
HTS dc SQUIDs (Sec.tV), and we compare those results 
with experimental data on HTS dc SQUIDs (Sec.|v|). 

II. ASYMMETRIC DC SQUID 

The asymmetric dc SQUID shown in Fig.|l| consists of a 
superconducting loop of inductance L intersected by two 
Josephson junctions with average values of critical current 
Iq, resistance R and capacitance C. The asymmetry in 
the junction parameters is described via the asymmetry 
parameters a/, an and ac, wich are defined according to 



4,i = Io(l 

I ,2 = J (l 



Oil) 



R 1 =R/{l-a R ), d = C(l-a c ), 



R 2 =R/(l + a R ), C 2 = C(l + a c ), 

(1) 

where the subscripts 1,2 denote the parameters of the left 
and right junction, respectively. 

Throughout this paper we consider two different origins 
of the junction asymmetry, noting, that in real devices we 
may have a combination of both: (i) geometric asymmetry 
and (ii) intrinsic asymmetry. 

In the case of geometric asymmetry we assume identi- 




Fig. 1 

The asymmetric dc SQUID 
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cal values of critical current density jo = Io/A, resistance 
times area p = R x A, and specific capacitance C = C/A 
for both junctions. Here, A = w x I is the junction area 
with width w and length I. We then introduce an asym- 
metry via different values for w, assuming constant I. In 
the case of bicrystal grain boundary junctions I equals the 
film thickness d, which can be assumed to be constant 
in practical devices. The geometric asymmetry is then 
described by the asymmetry parameter a g according to 
u>i = (l — a g )w and W2 = (l + a g )w, with w = {w\+w<z)/2, 
and we find the simple relation for the asymmetry param- 
eters a g — ai — an = a>c- 

In the case of intrinsic asymmetry we assume a g = 
and different values of jo and p for the two junctions which 
may reflect the natural spread in junction parameters. For 
simplicity we neglect the spread in C which is expected to 
be much smaller than the spread in jo and p. The intrinsic 
asymmetry is then described by the asymmetry parameters 
a.j and OL p according to 

10.1 =jo(l-ay), px=p/(l-a p ), 

V / 

10.2 = io(l + atj), p 2 = p/(l + a p ), 
and we get aj = a.j, an = a p , and ac ~ 0. 

1/2 

If we assume the scaling relation IqR = jop oc j as 
derived from the intrinsically shunted junction model |^| 
we derive the relations 

0L P = [1 - (1 - a?) 1 / 2 ]/^- or ctj = 2a p /{l + a 2 ). (3) 

Finally, if we consider a combination of geometric and in- 
trinsic asymmetry we can derive from the definitions of the 
asymmetry parameters given above the following relations 
between the asymmetry parameters 

a g ~ ai , a g -a a , . 

a J a/ttg-1' dI1U a P a R a g -l ' \V 

III. CRITICAL CURRENT VS. FLUX AND VOLTAGE VS. 
FLUX CHARACTERISTICS 

We derive now a very simple expression for ai, which 
can be used to determine its value experimentally without 
cutting the SQUID loop . Let us first consider the / c (^)- 
characteristics of the dc SQUID [c.f. Fig.0 for r = 0: the 
maximum critical I™ ax = / s + I 2 — 2Iq is maintained 
when I\ = 7o,i an d 1% = ^o,2- In this case the circulating 
current J is given as J(/™ aa: ) = (Jo, 2 — ^o,i)/2 = ailo, 
which is flowing in the SQUID loop if we apply an exter- 
nal flux $ + = oliIqL. Hence, the maxima of the I c (&)- 
characteristics are shifted by $ + along the $-axis, as com- 
pared to the symmetric SQUID with ai = 0, where I c is 
maximum at $ = 0. The same argument leads to a nega- 
tive shift $~ = —ailaL if the current direction is reversed. 
Hence, the maxima of the 7 C (^-characteristics for oppo- 
site polarity of the current arc shifted by A$ = $ + — $~ = 
2aiIoL. Using the screening parameter /3t = 2I L/$ we 
arrive at the very simple relation 
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Fig. 2 

Measured V($)-curves for different values of bias current 
I (positive and negative) at (a) T = 77K and (b) T = 63. 7K. 



which is also valid for finite values of T and for the shift of 
the minima of the U(<i>)-characteristics measured at con- 
stant bias current / w 2Io and —I rj — 2Iq. 

The clear-cut experimental determination of A$ requires 
the measurement of V($)-curves at various values of (3l, 
which can be varied by temperature. An example of such 
a measurement on one device at T — 77K and 63. 7K is 
shown in Fig.^. Plotting the measured flux shift, normal- 
ized by a/$o vs. @l should give according to Eq.(||) a linear 
dependence with slope 1, with the reasonable assumption 
that cei does not depend on T. Fig.|| shows the results 
of such measurements obtained for 6 different dc SQUIDs, 
where ai was obtained as a fitting parameter to give the 
expected slope of 1. For comparison, the results from sim- 
ulated U(<I>)-curves for various values of aj are also shown, 
which are in excellent agreement with Eq.(|^). 

Except for one device, these SQUIDs have intentionally 
been fabricated with a geometric asymmetry (a g ^ 0). 
From the known value of a g and the measured value of 
ai the asymmetry parameter otj can be calculated using 
Eq.([|). The results are listed in Table I. As a main result, 
we see that 3 SQUIDs show only a small asymmetry in jo 
with |afj | < 0.1. However, for the three other devices the 
asymmetry in jo is significant with values of |ay| up to 0.4, 
which demonstrates that the difference in critical current 
density for the two junctions can be quite large. 
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Measured normalized flux shift vs. /3_l for various dc 
SQUIDs. 



A$/$ = a//3 L , 



(5) 
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IV. TRANSFER FUNCTION: NUMERICAL SIMULATIONS 

As already evident from Fig. ^] the asymmetry can induce 
distortions of the V($)-curves, which leads to different val- 
ues of V$ and for the maximum positive and negative 
slope of the U($)-curves, respectively. To understand the 
impact of the asymmetry on the transfer function we per- 
formed numerical simulations to solve the equations for the 
phase differences 6x(t) and 62 (t) of the two junctions 
H over a wide range of values for /3l and the noise param- 
eter r = 27rfc j BT/Jo$o for both, the geometric and intrinsic 
asymmetry. In the latter case we assume the correlation 
between otj and a p as given in Eq.(|3|). 

Figure ^ shows V$ and V$ vs. 0l obtained for fixed 
T0l = L/L t h — 0.2 for geometric and intrinsic asymme- 
try. For ai = we closely reproduced the results ob- 
tained in fl for symmetric dc SQUIDs. In most cases 
we find for /3l~5 for the asymmetric SQUID a reduction 
of V$ as compared to the symmetric SQUID , which in- 
creases with decreasing 0l and increasing a/. However, 
in the case of geometric asymmetry, we find for interme- 
diate values of 1~/3l~5 an increase in V$ for ai=aa~0.5, 

which is concomitant with a strong distortion of the V($)- 

1/2 

characteristics. If we assume IqR oc j , we find that for 
geometric asymmetry an is always larger than for intrinsic 
asymmetry (for given a/). This implies that the distortion 
in V($) is dominated by the asymmetry in the junction 
resistances which becomes important for /3l~1. At 0^-0.2 
the reduction of and is similar for both types of 
asymmetry, indicating that for small 0l the asymmetry in 
the critical currents gives the main contribution to V$. 

Results similar to those shown in Fig]!] have been ob- 
tained over a wide range 1/80 < T0l < 0.5 which corre- 
sponds to 4pH< L <160pH for the SQUID inductance at 
T=77K where L t h = 321pH. We note that for the sym- 
metric dc SQUID it was shown in |jj that the normal- 
ized transfer function = V$$o/ IqR factorizes in — 
g{T0 L )-f{0 L ) where f(0 L ) is given as v^{0 L ; T0 L = 1/80). 
As a main result of our simulations for the asymmetric 
dc SQUID we find a similar factorization, with f{ai,0£) 
being now also dependent on a/, while g{Y0L) shows no 
dependence on a/. For this is shown in Fig.||(a) for 
geometric and in Fig.||(b) for intrinsic asymmetry. The 
simulation data shown in Fig.[5| can be approximated as 



(7.18 



8.42 

0.31 Q ' 



0.45 c 



1 \ r>0-lS+0A9aj' 



gtntAm) = [(8or / 3 i )°- 4 + o.9i(4r/? i ) 1 - 83 ]- 1 

(7.18-0.23^ 1 - 53 .0.4W- 7 )/3°- 18+0 - 46tt? " 



O-+0L) 



(6) 

which is plotted in Fig.[| as dotted lines. These equations 
enable one to calculate V£ immediately for any value of 0l, 
T(3l and a/, however only within the range of parameters 
displayed in Fig.||. Similar behavior was found for V { 
again with the factorization vZ 



<i> ■ 



g-(Tp L ) ■ f-(ai,0 L ) 



we did not derive expressions as 
As discussed above for /~ the 



with g~ w g + . For f~ 
for /+ given in Eqs([|) 
resistance asymmetry plays an important role. Its impact 
on V<$> needs to be studied in more detail. 

V. TRANSFER FUNCTION: NUMERICAL SIMULATION VS. 
EXPERIMENT 

To test the numerical simulation results we determined 
U* for various YBCO dc SQUIDs (#2-6 from Table I) with 
L w 2 and T ks 0.04 at T=77K. Figure | shows that the 
simulation (open symbols) predicts a reduction of both, 

and as compared to the symmetric case. For 
the experimental data (solid symbols) lie between the val- 
ues predicted for the symmetric and asymmetric case. Fol- 
ey/ = 0.4 the measured and calculated values for v~f are 
identical and lie almost a factor of two below the value pre- 
dicted for the symmetric SQUID. This supports the idea 
that at least part of the discrepancy in V$ between calcu- 
lated and measured values can be attributed to asymmetry. 



gteoPpL) 



[(80r/3 L ) a4 + 0.35(4r/3 L ) 



2.5i 




TABLE I 

Asymmetry parameters for various YBCO dc SQUIDs. 
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Fig. 4 

Calculated normalized transfer function v£ and V"^" vs. /3 l 

FOR r/?£, = 0.2 WITH GEOMETRIC ASYMMETRY (a), (b) AND 
INTRINSIC ASYMMETRY (c), (d); aj = (□), 0.2 (o), 0.4 (A), 0.6 
(V), 0.8 (O), 0.9 ( + ); SOLID LINES ARE GUIDE TO THE EYE. FOR 
COMPARISON. SIMULATION DATA (*) TOGETHER WITH THE ACCORDING 
FIT-FUNCTION (DOTTED LINE) FROM [1] FOR SYMMETRIC DC SQUIDS 
ARE ALSO SHOWN. 
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Fig. 5 

Calculated normalized transfer function v£ for (a) 

GEOMETRIC AND (b) INTRINSIC ASYMMETRY. THE PLOT SHOWS THE 
FUNCTION g{VPi,) WHICH DOES NOT DEPEND ON aj. THE INSET 
SHOWS THE NORMALIZATION FUNCTION /(a/,/?i). DASHED LINES ARE 
CALCULATED FROM EQ.([j). 



However, most measured values for are clearly larger 
than predicted for either the symmetric or the asymmetric 
SQUID. The reason for this deviation is most likely due to 
the large uncertainty in the value for an, which induces an 
increasing distortion with increasing an Q , Q . 

We finally note that a clear-cut comparison between sim- 
ulation and experimental data for V$ requires the exper- 
imental determination of at least a/, and preferably also 
of olr for a wide variety of HTS dc SQUIDs. Unfortu- 
nately, an extensive collection of such data does not exist 
yet. In order to obtain at least some information on the 
importance of asymmetry as a possible source for the dis- 
crepancy in V$ between theory und experiment we show 
in Fig.fi] the prediction for symmetric dc SQUIDs (dotted 
line) compared with experimental data taken from the lit- 
erature In addition, the large open symbols show our 
simulation results obtained for various values of ai , assum- 

1/2 

ing intrinsic asymmetry with IqR oc j . Obviously, the 
observed deviations between experiment and simulation for 
symmetric SQUIDs cannot be explained by asymmetry for 
large values of 0l~5. For smaller values of /?£, however, 
most of the experimental data lie within the range of val- 
ues covered by simulations which take into account asym- 
metry in the SQUID, although a large reduction of V$ due 



to asymmetry, say by a factor of five requires a very large 
ui k, 0.9. 

VI. Conclusions 

We have analyzed the performance of asymmetric HTS 
dc SQUIDs both experimentally and by numerical simu- 
lation, with focus on transfer function. Our simulations 
show that strong critical current asymmetry which may 
arise from a large spread in critical currents in HTS Joseph- 
son junctions can significantly reduce V$ for small (3l~2. 
This observation is important, since optimum performance 
requires the realization of small fit, w 1. We wish to stress 
that the asymmetry, which is most likely present in almost 
all HTS dc SQUIDs, may be one source for the previously 
found discrepancy in V$ between experiments and simula- 
tions, however it is not likely that this asymmetry is the 
major source of this discrepancy. 
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Comparison of measured and calculated transfer function. 
The asymmetry parameter ai was determined from the shift 
of the V($)-characteristics. For comparison the calculated 
values of v$> for symmetric dc squids are also shown. 
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Comparison of calculated and measured normalized 
transfer function plotted vs. f3 l . 
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